Acoustic fluidization for earthquakes? 
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Abstract : Melosh [1996] has suggested that acoustic fluidization could provide an 
alternative to theories that are invoked as explanations for why some faults appear 
to be weak. We show that there is a subtle but profound inconsistency in the theory 
that unfortunately invalidates the results. We propose possible remedies but must 
acknowledge that the relevance of acoustic fluidization remains an open question. 



1 Introduction 



In the standard rebound theory of earthquakes, deformation elastic energy is progres- 
sively stored in the crust and is suddenly released in an earthquake when a threshold 
is reached. The Ruina-Dieterich friction laws [Dieterich, 1972; 1978; Ruina, 1983] 
constitute the basic ingredient used to describe the interaction between the two sides 
of the sliding fault. Friction coefficients based on laboratory experiments [Byerlee, 
1977; Scholz, 1998] fail to account for modern observations of strain [Jackson et ai, 
1997], stress [Zoback et ai, 1987; Zoback, 1992a; 1992b] and heat flows [Henyey and 
Wasserburg, 1971; Lachenbruch and Sass, 1980; 1988; 1992; Lachenbruch et ai, 1995] 
(see Sornette [1999] for a synthesis). 

Resolutions of these paradoxes usually call for additional mechanisms, involving 
fluids [National Research Council, 1990], crack opening modes of slip [Brune et ai, 
1993], dynamical collision effects [Lomnitz, 1991; Pisarenko and Mora, 1994], fric- 
tional properties of a granular gouge model under large slip [Scott, 1996], space fiUing 
bearings with compatible kinematic rotations [Herrmann et ai, 1990], hierarchical 
scaling [Schmittbuhl et ai, 1996], etc. 

Melosh [1996] has recently suggested that the mechanism of "acoustic fiuidiza- 
tion" could provide an alternative to theories that invoke pressurized fluids as an 
explanation for why some faults appear to be weak. Fluidization usually refers to the 
experimental observation that granular material in the presence of an interstitial fluid 
can liquidify when shaken sufficiently strongly [Russo et ai, 1995]. The liquifaction is 
due to the fact that granular media are first compressive for small deformation lead- 
ing to an increase of the interstitial fiuid pressure. This increase in turn decreases the 
friction between the grains that can eventually become free to shear. For example, 
liquifaction of sediments by resonant amplified seismic waves have been proposed to 
be in part responsible for the damage and collapse of certain buildings during the Mi- 
choacan earthquake, 1985 [Lomnitz, 1987] and for the damage in the Marina district 
of San Francisco during the Loma Prieta earthquake [Bardet, 1990]. 

In the acoustic fluidization mechanism [Melosh, 1979; 1996], no interstitial fluid is 
invoked. A fraction e of the earthquake energy is released as high-frequency acoustic 
waves that scatter off and shake the granular gouge leading to the build-up of a 



local acoustic pressure. According to [Melosh, 1996], when this pressure becomes 
of the order of the overburden lithostatic pressure pgh, the granular gouge becomes 
essentially free to slip without much residual friction. 

The purpose of this note is to show that there is a problem with this mechanism 
because it predicts a slip velocity during an earthquake more than two orders of 
magnitudes smaller than the typical meters per second for observed earthquakes, in 
contradiction to the result of Melosh [1996]. The problem stems from a confusion 
in the definition of dissipation and scattering lengths. We then suggest possible 
modifications of Melosh 's theory that could resolve this problem and which lead to a 
richer theory. 

2 Summary of Melosh's theory and useful back- 
ground 

2.1 Acoustic wave energy transport 

The first ingredient is the generation and transport of high-frequency acoustic waves 
in the core of the fault. Melosh [1996] uses the standard diffusion equation (his 
equation 2) for the elastic transport of acoustic waves [Ishimaru, 1978; Sornette, 
1989a-c] with a dissipation and a source term : 

where E is the acoustic wave elastic energy density. The diffusion coefficient ^/4, 
where { is called the scattering diffusivity by Melosh [1996], can be expressed in 
terms of the elastic mean free path 1^ and of the transport acoustic wave velocity c 
at scales below Ig. The velocity c is of the order of the shear wave velocity [Turner 
and Weaver, 1996; Van Albada et al, 1991; Van Tiggelen and Lagendijk, 1993]. We 
thus get [Ishimaru, 1978; Sornette, 1989a,c] : 

l-ici.. (2) 

We stress that the term "elastic" refers to the fact that le is the characteristic distance 
over which an acoustic wave propagates before being scattered in other directions, 
without any loss of energy. 



The l.h.s. and first term of tfie r.fi.s. of (|lD give the diffusion equation which 
describes the transport of an acoustic wave in a multiple-scattering medium. The 
second term of the r.h.s. of (|l]) will be shown to describe the presence of a genuine 
absorption, while the last source term corresponds to the conversion to acoustic waves 
of a fraction e of the mechanical work performed per unit time by the fault motion 
with strain rate e and shear stress r. 

2.2 Diffusive transport 

The first two terms of 

dE/dt = {^/4)V^E (3) 

gives the standard parabolic diffusion equation which is based on the following pro- 
cesses. Once generated from a source, an acoustic wave propagates roughly ballisti- 
cally over a typical distance of the order of the elastic (scattering) mean free path 
le- Over this distance, the equation governing the acoustic wave propagation is the 
hyperbolic wave equation 

where the wave amplitude A is related to E hj E = |Ap and the wave velocity c may 
depend locally on position to reflect the heterogeneity of the medium. Due to this 
heterogeneity, the wave is scattered off its initial propagation path along the direction 
X and its intensity in this direction x decays as exp[— x/Zg]- This exponential decay of 
the intensity does not correspond to a genuine absorption but rather reflects the loss 
of acoustic energy along the direction x to all possible scattered waves in all other 
directions. Mathematically, the exponential decay exp[— x//e] can be derived from 
(§) using standard scattering theory [Ishimaru, 1978]. The conservation of acoustic 
energy is ensured by the fact that the sum of wave intensity over all directions of 
propagation remains constant. 

Beyond the distance le, the nature of the transport of the wave intensity crosses 
over from ballistic (i.e. straight propagation) to diffusive, corresponding to the picture 
where the acoustic wave can be viewed as a superposition of random walks with 
typical step length equal to le- This means that the equation for the wave propagation 
changes from the hyperbolic wave equation for the wave amplitude to the parabolic 
diffusion equation for the wave intensity given by the first two terms of 



One can quantify this by the following example. Consider an acoustic wave of 
energy Eq impinging on a slab of thickness L made of heterogeneities that scatter off 
the acoustic wave, and whose scattering strength is quantified by the elastic mean free 
path /g- Anderson [1985] and Sornette [1989c, d] have revisited this diffusion equation 
to get the transmission coefficient in this example, i.e. the acoustic energy which is 
transmitted to the other size of the slab, as a function of its thickness L. The result 
is : 

EiL) - I • (5) 

Note that the decay follows the algebraic 1/L law rather than an exponential law. 
Furthermore, the acoustic intensity profile within the slab is linear and not exponen- 
tially decreasing : 

E{z) ^ Eo ^^y ~^ , ioY l,<x<L-k. (6) 
ij 

These results (HJ^) highlight that the diffusive transport of the acoustic energy due 
to multiple scattering event is very different from the exponential attenuation that a 
genuine absorption would produce. 

2.3 Absorption 

The third term —{c/\Q)E of (0) quantifies genuine absorption processes. The pa- 
rameter A is the acoustic wavelength and Q is the quality factor. To see that this 
term reflects absorption, we consider (|l[) in absence of the spatial derivative V^-E and 
of the last source term : 

— = -—E . (7) 
dt \Q ^ ' 



Its solution is 



E(t) = E„ exp 1-^ t 



which is very different from the energy decay (^) solely due to diffusion. It is thus 
clear that the term —{c/XQ)E is not coming from elastic scattering but solely from 
genuine absorption, i.e. conversion of acoustic energy into thermal energy. 
The usual definition of the quality factor Q is [Knopoff, 1964] 

Q ^ 2vr^ , (9) 



where la is the absorption length defined by the exponential decay exp {—x/la) of a 
hallistically propagating wave in an absorbing medium. 

Melosh [1996] introduces a characteristic length which he calls (misleadingly) 
the "scattering length" , defined by 



Using @) and (||), we get 



\flla. (11) 

Calling a "scattering length", is misleading because is in reality the effective 
absorption length in the diffusive medium. To see this, we use the standard diffusion 
relation 

ll ^ 6|r, (12) 

linking the radius of gyration covered by a diffusing process over a time Tq = la/c 
equal to the time needed for the wave to cover the real distance along its convoluted 
multi-scattered path. The prefactor 6 holds for diffusion in a three dimensional space. 



Using (^, we get ^ v^2lj^, which recovers the (|Tl]) up to a numerical factor of 
order unity. The expression (|TT]) can be derived by several other methods [Sornette, 
1989d]. What is important is that scales as the geometrical mean of le and la, 
which comes from the random walk nature of the diffusive process. Physically, in the 
diffusive regime, the acoustic wave energy is absorbed over the characteristic length 

which stems from the fact that, to cross the distance the wave follows random 
walk paths of length la ~ ll/le- This reflects that attenuation of a wave in a scattering 
medium is a function of both absorption of energy and scattering. 

2.4 Feedback of the acoustic vibrations on the sUp rate 

The interesting idea of Melosh [1996] is that the high-frequency vibrations may shake 
the fault and unlock it, leading to an easier sliding motion. For this, he proposes the 
following effective friction equation, relating the strain rate e, the shear stress r and 
the normalized acoustic wave energy = E pc^ /{pghY '■ 

. r r i-^^f(i7i) ^ 

pAc [l + erf(^). 



(13) 



The main physical phenomenon taken into account in this equation is that, due to 
the acoustic shaking, the effective viscous friction r/e is a decreasing function of the 
acoustic wave energy. This mechanism is related to the velocity weakening mechanism 
induced by collision between asperities that lead to a transfer of momentum from the 
direction parallel to the fault to the direction transverse to it [Lomnitz-Adler, 1991; 
Maveyraud et al, 1998]. 

Putting ([T3|) in ([l|) and looking for stationary modes gives the non-linear ordinary 
differential equation (^) (for the case ?7 = 1, see below), whose analysis leads to the 
prediction of two rupture modes [Melosh, 1996]. 

3 Problem with Melosh's theory 

In order to obtain realistic values, there are some constraints that the model param- 
eters must satisfy. The key parameter is the "regeneration" parameter 



where S = r/pgh is the normalized shear stress. Melosh [1996] finds reasonable 
solutions only for rS^ > 2.8. For a typical fraction e ~ 0.1 of conversion to acoustic 
waves of the mechanical work performed per unit time by the fault motion and for a 
ratio r/pgh as low as 0.1 as suggested from observations on the San Andreas fault, 
this value rS^ > 2.8 corresponds to Q > 5600. This estimation may vary by an order 
of magnitude with the conversion factor and the relative shear stress. However, the 
message is that the quality factor Q measuring the attenuation of the acoustic waves 
must be high, in the range of 10^ for the acoustic waves to be self-sustained during 
the earthquake slip motion. This is the first condition. 

On the other hand, Melosh's theory predict the slip velocity 



during a typical earthquake. Using a shear stress r ^ 10 MPa, cl density p = 
3000 kg/m? and c = 4 km/ s gives u ~ 1.2 (/*/A) m/s. Thus, a realistic slip velocity 
u ^ 1 m/s requires that 




(14) 
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(15) 



(16) 



Together with (|TTD and (|), this leads to 

Q ^ 3 ^ . (17) 

'e 

This last expression (|T^ is totally incompatible with the above condition Q > 10^, 
as this would lead to /e ~ A/100 or smaller. This last condition is a physical impos- 
sibility: the elastic scattering length is always much larger than or at the extreme 
limit of the same order as the wavelength. The physical intuition is that a wave is 
defined over a length scale of the order of the wavelength (otherwise, there are no 
spatial oscillations) and the scattering process needs at least this scale to operate. 
The limit /g — A is attained only under exceptional circumstances leading to a novel 
phenomenon, called Anderson localization, in which the acoustic wave do not prop- 
agate anymore but oscillate locally. Extraordinary efficient scatterers are needed to 
reach this regime [Sornette, 1989c]. It is thus clear that the condition Zg ^ A/100 is 
utterly unphysical. 

If in constrast, we put Q ^ 10^ in ([TOD , we get l^, ~ 160 A, which from ([TSD leads 
to a maximum slip velocity ii ~ 7.5 mm/ s, using the numerical example of Melosh 
[1996]. This slow sliding velocity is unrealistic for earthquakes. 



4 Possible remedies 

A first remedy is to relax the condition used by Melosh that the acoustic pressure 
needs to reach the overburden pressure in order to significantly affect the fault friction. 
We propose that only a small fraction rj of it is enough to liquidify the fault. 

Indeed, it is well-established experimentally [Biarez and Hicher, 1994] that the 
elastic modulii of granular media under large cyclic deformations are much lower than 
their static values. This effect occurs only for sufficiently large amplitudes of the cyclic 
deformation, typically for strains above 10~^. At = 10~^, the elastic modulii 
are halved and at ea = 10~^, the elastic modulii are more than five times smaller 
than their static values. As a consequence, the strength of the granular medium is 
decreased in proportion. Melosh (private communication) also finds in laboratory 
experiments that a large decrease in elastic modulus is required to fit the measured 
fiow rate of acoustically fiuidized debris. This is consistent with fiow in granular 



material in a completely dilatent state and agrees with measurements reported in the 
literature for the elastic modulus of highly strained granular materials. 

In absence of cohesion forces, the strength of a granular material is solely due to 
the effect of gravity weight that put grains in contact together and the resistance to 
shear is governed by Coulomb's law according to which the shear stress at the thresh- 
old for sliding is proportional to the normal stress. The coefficient of proportionality 
is the friction coefficient. As we have discussed above, in presence of acoustic flu- 
idization, the resistance to shear deformation is decreased as a consequence of the 
reduction of the effective elastic modulus. Correlatively, the threshold for sliding 
is also decreased (strength is often proportional to elastic modulus in brittle mate- 
rial). In order to capture these phenomena, we propose, what is maybe the simplest 
approximation, that the criterion for unlocking the fault is changed from Melosh's 
criterion to the condition that the acoustic pressure needs only reach a fraction t] 
of the overburden pressure. Let us stress that the essence of our argument is that 
the acoustic energy is fed by the moving fault and thus the acoutic particle velocity 
adjust to the changing elastic modulus so that the overall acoustic energy is "exter- 
nally" controlled by the rate of global elastic dissipation. The acoustic fluidization 
thus controls the shding threshold rather than solely the acoustic particle velocity. 

We need to estimate the strain created by the acoustic field. The acoustic pressure 
is related to the acoustic particle velocity v by 

p — pcv . (18) 

Assuming 

P = vpgh , (19) 

this yields 

V = righ/c ^ 12 m/s (20) 

for p fa 200 MPa, a density p = 3 10^ kg/m^, a velocity c = 4000 m/s and r] — 0.1. 
This corresponds to an acoustic wave displacement Ua — f/27r/ fs 2 10~^ m at a 
frequency / ~ 10^ Hz. The corresponding strain Ua/w is ~ 2 10~^ for a gouge width 
w of the order of one meter [Melosh, 1996] over which the intense shaking occurs. 
These estimations suggest that Melosh's criterion that the acoustic stress fluctuations 



must approach the overburden stress on the fault for acoustic fluidization to occur is 
too drastic and smaller shaking can reduce significantly the fault friction. 

Persuing this reasoning, we see that the fundamental equation (6) in [Melosh, 
1996] is changed into 



where rj = 1 recovers the case treated by Melosh. \E' is the normalized acoustic energy, 
( = z/l^, z is the coordinate perpendicular to the fault, the regeneration parameter is 
r = eQ/2 where e is the acoustic energy conversion efficiency, S = r/pgh and erf(x) 
is the error function. We see that a factor t] < 1 implies a more effective generation 
of acoustic waves because the second source term of the r.h.s. is larger. But, since 
the bracket term saturates to one for large energies and/or small 1], this does not lead 
to a significantly larger slip velocity than found above. This remains a problem of 
the theory. 

This problem might be alleviated by treating e self-consistently as a decreasing 
function of the friction coefficient, and thus of the acoustic energy density. The 
problem then becomes even more non-linear because it reflects in addition the de- 
pendence of the acoustic radiation efficiency of the granular gouge on the amplitude 
of the acoustic vibrations. In addition, the elastic modulus is also really nonlinear 
and it is only the tangent modulus that decreases close to the yield in the dilatent 
region, which suggests that the above linear formalism is not adequate and should be 
modified. 

Further improvement could also take into account that the stochastic acoustic 
energy may deviate from a Rayleigh distribution [Ishimaru, 1978; Mirlin et ai, 1998] 
when the medium is strong heterogeneous. This modifies the functional form of the 
term in bracket in eq. (|2l|) and thus all numerical estimations. 

We ackowledge stimulating discussions and correspondence with H.J. Melosh. 
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